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We review a performance of Fock space methods in calculating spectra of a range 
of supersymmetric models with gauge symmetry. Examples include: a) SU(2) Su- 
persymmetric Yang Mills Quantum Mechanics in four euclidean dimensions, b) 
Quantum Mechanics of one fermion and one boson with infinite number of colours, 
and c) planar 1+1 dimensional Yang Mills theories with adjoint matter. Infrared 
divergencies of the latter theories with scalars are briefly discussed and a possible 
dynamical solution of the problem is suggested. 



1 Quantum mechanics and a space reduced field theory 

It is usually believed that simple quantum mechanics of the finite number of degrees 
of freedom can hardly teach us something about the quantum field theory with all its 
subtleties like renormalization, spontaneous symmetry breaking, etc. This is in general 
correct, however there exists a class of quantum mechanics which are being intensively 
studied precisely because of their connection to the field theoretical systems [TJ [21 13]- 
Consider the following hamiltonian i = 1, ..,D — 1, a = 1, TV 2 — 1. 

H = \^ a + 9 ^f ab cfa d eXlx^ d xi + ^f abc ^T k ^ b xt (1) 

It describes a quantum mechanical system with the finite number of degrees of freedom 
(e.g. 15 for D=4 and N=2) and results from the dimensional reduction of the D dimen- 
sional supersymmetric Yang-Mills theory to one point in the D-l dimensional space. In 
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spite of the reduction, it still has many nontrivial properties inherited from the parent, 
space extended, theory [4J, [5l E] and can teach us quite a lot about the latter. Even the 
SU(2) model is not soluble for D=4. We have studied it numerically by constructing a 
gauge invariant basis of bosonic and fcrmionic Fock states, cutting off the total number 
of bosonic quanta and subsequently increasing the cutoff [7] . The algorithm worked very 
satisfactorily and we were able to uncover a rich, manifestly supersymmetric spectrum 
with coexisting localized and non- localized states, SUSY vacua and a fractional bulk 
value of the Witten index [5J [TO] . All these results agree with [TT] , and in some cases 
extend [3], theoretical predictions. 



2 Large N quantum mechanics 

There is much interest in the large N limit of quantum systems and above approach turns 
out to be very useful in studying such models as well. Direct calculations allow to obtain 
spectra for the first few lowest values of N. An extrapolation to N = oo does not seem 
feasible, however. Nevertheless, it appears that one can calculate analytically matrix 
elements of typical hamiltonians directly at N = oo. |12J . We have studied with this 
method a close cousin of the space reduced SYM^ in the planar limit. Its hamiltonian 



reads H = {Q, Q*}, Q = y/2Tr[fat(l + gat)],Qi = V2Tr[p(l + ga)a], or explicitly 



It was found that, in spite of its simplicity, the system exhibits many interesting phenom- 
ena: unbroken supersymmetry, the phase transition in the 't Hooft coupling, at A = 1, 
exact duality between the strong and weak coupling phases, rearrangement of supermul- 
tiplets across the transition point and emergence if the new vacua in the strong coupling 
phase [HI EH E] • All this was first discovered numerically and subsequently confirmed 
by the analytic solution. Moreover, at strong coupling the model is exactly equivalent 
to the XXZ Heisenberg spin chain and, at the same time, to the lattice gas of q-bosons. 
This proves existence of a hidden supersymmetry in these well explored statistical models 
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Tr[a^a + g(a* a + a f a 2 ) + g 2 a t2 a 2 ] 
Tr[pf + g{Pf{(J + a) + /V + o)f) 
5 2 (/t a /at + /W/ + /t/ a t a + /t a t/a)]. 



nans]. 
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3 Planar field theories in 1+1 dimensions 

Extension to the field theoretical systems is in principle straightforward. One has to 
diagonalize the Hamiltonian matrix calculated in the physical basis of harmonic oscilla- 
tors. The obviously crucial difference is that now we deal with the infinite number of (e.g. 
momentum) degrees of freedom. Nevertheless one can define cutoff schemes which allow 
to extract meaningful continuum physics. There exist two popular ways to introduce a 
cutoff. 

Light Cone Discretization (LCD) [T7] replaces the total momentum P of a proton, 
say, by an integer K. This momentum can then be split between various numbers of 
partons, each carrying, an integer momentum r% > 0. Therefore all partitions of K 
into sets of integers {r} define a finite Fock space bound by only one cutoff, K, which 
discretizes momenta and, at the same time, cuts the multiplicities of partons [17 ] [18 | IT9] . 
The Hamiltonian matrix ({r}|ii"|{s}) is then calculated and diagonalized, similarly as in 
the case of quantum mechanics. The whole art consists of the meaningful extrapolation 
with K — > oo. 

The second approach consists of solving the Integral Equations (IE) [20] which are 
equivalent to the eigenequation H\&) = M 2 |$) (in the (LC) formulation a Hamiltonian 
is proportional to a mass 2 operator). Decomposing a bound state into its Fock com- 
ponents |$) = X^2/[^ x ]^(l — X1 — X2 — .. .x„)<& n (xi,X2, ■ ■ .x n )\xi, . . . , x n ) , turns the 
eigenequation into an infinite hierarchy of integral equations 

) = A(8)$„ + S®$„_ 2 + C*(8)$„ + 2, (2) 

where each term is a convolution of the wave function with the amplitude for scattering, 
emission and fusion of partons respectively. Appropriate amplitudes can be readily read 
from the LC form of the hamiltonian. In practice one has to cut the hierarchy limiting 
the multiplicity of partons. At fixed maximal multiplicity n max the LCD is indeed 
the momentum-discretized version of IE cut to n < n max . In Fig.l we compare the 
LCD simulations restricted to n nmax = 2 with the IE solution at the same multiplicity. 
The agreement is satisfactory, however very fine discretization is required to see the 
convergence (K ~ 2000)0. This is due to the singular nature of the scattering amplitude 
( A ~ Pjj ). If one allows for arbitrary number of partons in the LCD approach the 
rapid growth of the number of states with K excludes K > 25 which makes extrapolation 
to K = 00 rather delicate (151(151120] . 

In contrast solving integral equations with two partons required only few basis functions 
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Figure 1: Comparison of the LCD and IE results (symbols at 1/K = 0) 



4 Infrared divergencies for theories with scalars 

The previous discussion applies to a generic LC Hamiltonians. The calculations reported 
in Fig. [3] were done for the YM2 with adjoined fermions. There is also some interest in 
YM2 with adjoined scalar matter (which can be thought of as the dimensional reduction 
of YM3), and finally one may combine both into the two-dimensional supersymmetric 
Yang-Mills theory, SYM2 |21j . However introducing scalars brings one more complica- 
tion: the integral equations have infrared logarithmic divergence. It appears explicitly in 
the mass term which is part of the diagonal (in multiplicity) transition in . One way to 
dispose of it is the " mass renormalization" introduced in [20j . This might be possible for 
the YM2 with a scalar matter, however for the supersymmetric model such a procedure 
would break supersymmetry. The attractive possibility to maintain SUSY would occur 
if the above divergence was canceled dynamically by other contributions [22] . 

Some support for this idea is provided by Fig. 2 where we compare the lowest mass 
obtained by the LCD simulations with (open circles) and without (filled circles) the 
multiplicity cutoff. Clearly the dependence on the cutoff is weaker (indicating possible 
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Figure 2: The LCD results without (filled) and with the multiplicity cut (open,n r , 
2). 



convergence at larger K) when all multiplicities arc allowed. Second piece of evidence 
comes from careful inspection of the integral equations @ which shows that, when the 
wave functions have a suitable divergence at x ~ 0, additional divergences occur which 
in fact cancel the original IR divergence of the mass term [22) . For this mechanism to 
work one has to include all Fock components in eqns ©■ This is consistent with the 
message learnt from Fig. 2. In fact the whole mechanism is analogous to the classical 
Bloch-Nordsieck treatment of the IR singularities in QED as discussed in the following 
Section. 
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5 The Bloch-Nordsieck inspired toy model 

A good insight into the mechanism of the above dynamical cancelations is provided by 
the following Light Cone model hamiltonian 



f P 9 
H = / dk [atfcOfc - jk(a\ + a k ) + j%\ = V" H k . j k = — = 

Jo k Vk 



(3) 



H k =A\A k , A k = e- iPkjk a k e iPkjh , P k = -±=(a k - a) k ). (4) 

?V 2 

The model is of course soluble, the eigenstates are the BN coherent states 

\n k ) new = —f=\0)new = er %Pk3k \n k ) oW (5) 

Vn! 



and the eigenvalues are integer M 2 = J2 k n k- 

The integral equations ((2|) can be easily derived 



1 

2, 



M f n (x 1 , . ..,x n ) = I n + y ^ ) • 

<_1 „ ^0 



(6) 



and are divergent for the sources ([3]), but in fact they must give (and they do!) the above 
finite spectrum. 

We expect that the similar cancelations occur in the case of YM2 with adjoint scalars 
and also in SYMi- Figure 3 lends some support for this analogy On the left hand we 
show numerical solutions of ([6]) for the first three eigenvalues as the function of the IR 
cutoff e. Three curves for each eigenvalue correspond to increasing multiplicities (top to 
bottom: n max — 2, 3, 4). At finite n max eigenvalues are divergent at small e as expected. 
However for increasing n max the divergence is shifted towards smaller and smaller e 
and the integer values of M 2 are better and better approximated. On the right hand we 
compare in the same way the masses obtained from the LCD and IE for YM2 with adjoint 
scalars 22 . Both calculations were done only in the two-parton sector. Results are not 
far from each other showing the onset of the divergence at low e. It is expected that as 
we include higher Fock sectors the two curves would converge and the IR divergence will 
shift to yet smaller e, similarly to the toy model example. 
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In the summary, reliable extrapolation of LCD data requires very large cutoffs, while 
the IE do not have this problem. On the other hand, LCD samples better multiparton 
Fock states which is difficult to achieve with the IE. In that sense the two methods are 
complementary and only comparison of both can provide unbiased information about 
many-parton phenomena in the continuum limit of quantum field theories including QCD. 
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